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ABSTRACT 



J>" We present explicit form of non - adiabatic coupling (NAC) elements of nuclear Schroedinger equa- 

^ tion (SE) for a coupled tri - state electronic manifold in terms of mixing angles of real electronic 

O basis functions. If the adiabatic - diabatic transformation (ADT) angles are the mixing angles of 

"nJ" electronic basis, ADT matrix transforms NAC terms to exactly zeros and brings diabatic form of 
O 

""^ SE. ADT and NAC matrices satisfy a curl condition and find a novel relation among the mixing 
• rH angles for irrotational case. We also find that extended Born - Oppenheimer (EBO) equations have 

^meaningful solution and can reproduce numerically exact results only when the equations are gauge 

P-invariant. 
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Herzberg and Longuet - Higgins (HLH) [1] made an interesting observation in molecular dynamics 

- a real valued electronic wavefunction changes sign when the nuclear coordinates encircles a closed 
path around a conical intersection. This is so - called geometric phase (GP) effect, where the total 
wavefunction is not well behaved and the Born - Oppenheimer (BO) treatment [2] breaks down. 
HLH corrected this problem in an ad hoc way by multiplying a complex phase factor with real 
double - valued electronic wavefunction such that the resulting complex electronic wavefunction is 
single valued. Mead and Truhlar [3] generalized the BO equation by introducing a vector potential 
into the nuclear Hamiltonian to account this effect and the approach is reminiscent of HLH complex 
phase factor treatment. The effect has been found to show up immense significance, particularly, 
in scattering processes [4]. 

The BO treatment is based on the fact that slow - moving nuclei can be distinguished from fast 

- moving electrons in molecular systems. This distinction leads to the BO approximation, which 
includes the assumption that NAC terms are negligibly small and the upper electronic states do 
not affect the wavefunction on the lower state. Though the implication of this assumption is 
considered to be independent on the energy of the system, the ordinary BO equation was also 
applied for cases with large NAC terms assuming that upper electronic states are classically closed 
at low energies. Even if the components of total wavefunction on the upper electronic state(s) are 
negligibly small at enough low energies, their product with large NAC terms are definite values and 
the BO approximation breaks down for those cases. 

Since NAC terms appear as off - diagonal elements in the SE [5], formulations of EBO equation 
are the attempts [6-7] to include the effect of off - diagonal (non - adiabatic) elements on the 
single surface equation. Baer [7] derived a new set of potential energy coupled BO equations for 
two electronic states, where the effects of NAC terms are being translated on the diagonal. At 
enough low energy, Baer assumed that upper electronic state is classically closed and approximate 
EBO equations for both the surfaces were formed. This EBO equation was used to calculate 
transition probabilities in a two - arrangement - channel model [8-9] and reproduced the correct 
results obtained from rigorously derived diabatic SE. Varandas and Xu [10] expressed NAC elements 
of adiabatic nuclear SE in terms of nuclear coordinate dependent mixing angle of two real electronic 
wavefunctions, found that ADT angle [11] is identical (up to an additive constant) with the mixing 
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angle and indeed, formulated single surface EBO equations in the vicinity of degeneracy. 

Baer et al. [12-13] considered coupled BO equations for N (> 2) adiabatic nuclear SE and derived 
single surface EBO equations for a model situation ( not a general case). These EBO equations 
contain the effect of NAC terms where the ground adiabatic PES have a degeneracy with excited 
surfaces along a single line. Even though this derivation was not persuaded for a general situation, 
the formulation shows that such an extension is viable and those EBO equations will have meaningful 
solution only when they are gauge invariant. 

In this article, we carry out the BO treatment of a coupled three state electronic manifold from first 
principles. The BO expansion of molecular wavefunction and the total electron - nuclei Hamiltonian 
in the adiabatic representation are given by: 



*(n,e) = ^^(n)&(e,n), 

i=i 

H = f n + H e (e,n), 
f) 2 

T = — — W 2 

H e (e, n)&(e, n) = u;(n)&(e,n), (1) 

where n and e are the sets of nuclear and electronic coordinates, respectively, £j(e, n) is the eigen- 
function of the electronic Hamiltonian, H e (e, n) with eigenvalue, ttj(n), also ^(n) is the nuclear 
coordinate dependent expansion coefficient and T n is the nuclear kinetic energy (KE) operator. 
When we substitute equation (1) in the time - independent SE and multiply by electronic ba- 
sis, £j(e, n), the matrix representation of adiabatic nuclear equation takes the following form after 
integrating over the electronic coordinates, e, 

3 

J2(Hij - £<yV» = 0, 2 = 1,2,3, 
J'=l 

Hu = -^(V 2 + 2ff -V + ^+^n), 



2m 

H * = -^(2^ ) -V + ^ ) ) = < 
Tf = <&(e,n)|V| (e,n)>, = <fc(e, n)| V 2 ^(e, n)), 

(^(e,n)|^(e,n))=^, (2) 



^ 2 ,, 
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We define the following three real orthonormal functions as the electronic basis set (£1, £ 2 and £ 3 ) 
where a(n), (3(n) and 7(11) are the nuclear coordinate dependent mixing angles (e.g. a(n) is the 
mixing angle between electronic states 1 and 2), 



^ — cos a sin (3 sin 7 — sin a cos 7 ^ 
— sin a sin (3 sin 7 + cos a cos 7 
cos /? sin 7 / 



^ cos ct cos (3 \ 
£1 = sin a cos /3 , £2 = 
\ sin /3 / 

' — cos a sin (3 cos 7 + sin a sin 7 \ 
£3 = — sin a sin f3 cos 7 — cos a sin 7 , 
\ cos [3 cos 7 / 

and rewrite the kinetically coupled nuclear SE (equation (2)) as below, 



(3) 



2m 
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where the NAC matrix (r(= r^)) is defined as, 

/ 



with matrix elements, 



T 



-t 





u 2 -E 
u 3 - E J 



^2 



r 3 



r 2 \ 



V -* 2 -* 3 ; 



0, 



(4) 



(5) 



ti = — cos (3 cos 7 Vet — sin 7V/3, 

t 2 = cos/?sin7Vo; — COS7V/5, 

_» — » — > 

i 3 = — sin /SVo; — V7. 



(6) 
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When we substitute * = A^ d in equation (4) with the following choice of ADT matrix, 



A = 



cos a cos (3 



sin a cos (3 



sin (3 



\ 



— sin a cos 7 cos a cos 7 cos j3 sin 7 

cos a sin /3 sin 7 — sin a sin (3 sin 7 



(7) 



sin a sin 7 — cos a sin 7 cos /? cos 7 

^ — cos a sin (3 cos 7 — sin a sin /3 cos 7 y 

adiabatic nuclear SE is being transformed to the potentially coupled diabatic nuclear SE, 

E{(-^V 2 -£)<% + ^}^ = 0, i = 1,2,3, 
W = A f UA, C/ij = li^. 



(8) 



Since we find that the above form of r (equation (5) and (6)) and A (equation (7)) matrices satisfy 
the ADT condition [11], 



VA + f A = 0, (9) 

we arrive an equation (known as curl condition) for each NAC element, 7^, considering the analyt- 
icity of the transformation matrix A for any two nuclear coordinates, p and q, 

—T g - - —t p - = (r q T p ) . . - (T P T q ) . . 

Tfi =< 6|v p |0 >, 4 =< 6|v 9 |0 > . (10) 

The curl condition for each pair of electronic basis, {|£i >, |£ 2 >}, {|£i >> I £3 >} an d {I62 >, (£3 >} 
is satisfied and the explicit forms of curl equations are the following, 



Curir 12 = [r x r] 12 = sin/3cos7[V p a:Vq/? — V q aV p (3] + sin7cos/3[V p o;V g 7 — V g aV P 7] 

- cos 7 [V p /?V ? 7- VqPVp-Y], 
Curir 13 = [r x r] i3 = — sin (3 sin 7[V p a;V g /3 — V q aV p (3] + cos 7 cos /3[V p aV g 7 — V g aVp7] 

+ sin7[V p/ 3V g7 - V,/3V P 7], 
Curlr 23 = [t x r] 23 = cos/?[V p a;V ? 7 - V,aV P 7]. (11) 
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Since r(n) goes to zero rapidly enough as the radial coordinate tends to infinity [14], r may be 
resolved into an irrotational and a solenoidal part [15-16]. On the otherhand, the explicit form of 
divr^s are given by: 



divf 12 = 2 sin/9cos/9sin 7(Va) 2 + 3 cos/3 sin7VaV7 — 3 COS7V/3V7 — cos/3cos7V 2 a 

— sin 7V 2 /? — sin/3cos7VaV/3 

divT 13 = 2sin/3cos/3cos7(Va) 2 + 3cos/?cos7VaV7 + 3sin7V/9V7 + cos/?sin7V 2 a 

— cos 7V 2 /? + sin (3 sin jVa V/3 

divT 2 3 = 2cos 2 /3sin7cos7(Va) 2 — 2sin7cos7(V/9) 2 — 3cos/9cos 2 7VaV/3 

+ cos ft sin 2 jVa V/9 - sin (3V 2 a - V 2 7 (12) 

Since Va, V/3 and V7 (generally V 2 a, V 2 /5 and V 2 7 also) are not zero in the vicinity of the 
conical intersection, divry 7^ for any value of mixing angles, i.e., vector field r corresponds to 
non - solenoidal case[16-17]. In the vicinity of the conical intersection, we presently handle only the 
irrotational part, i.e., curl equations are Abelian, CurWij = 0. Thus, equations (11) have unique 
solution as below (if (3 ^ f or ^ ^-), 



V p aV q f3 = V q aV p /3, 
V P /3V ?7 = V^V p7 , 

V p aVg7 = V 9 aV p 7, (13) 
with the implication that nuclear coordinate dependent mixing angles are related by integer ratios, 



a(n) : /3(n) : 7 (n) — k : I : m. (14) 

If an unitary transformation matrix, G (\1/ = G$), diagonalizes NAC matrix, r, with eigenvalues, 
and ±iu, the adiabatic SE (equation (4)) transforms as, 

h 2 - 

(V + iu) 2 <$> + (V - E)$ = 0, V = G*UG, 

2m 



io = ±^t 2 + t 2 + t 2 

= ±{(Va) 2 + (V/?) 2 + (V7) 2 + 2sin/5VaV7}i (15) 
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One can rewrite the product, V$, for the ith equation as, = + Y^=2^ij( u i ~ u i)' l Pj-, 

i — 1,2,3 and impose the BO approximation, j^l >> |^|, i = 2, 3 (considering that at enough 
low energy, both the upper electronic states are assumed to be classically closed) to form the single 
surface adiabatic nuclear SE [12], 

h 2 

- — (V + iuj i ) 2 ^ l + (u 1 -E)^ i = 0, i= 1,2,3. (16) 

Equation (14) simplifies both the adiabatic and diabatic equations (4) and (8), respectively for any 
arbitrary ratios of mixing angles, i.e., the NAC matrix takes the form, r = Va.g(a), where g(ct) is 
mixing angle dependent 3x3 matrix. At the same time, we know that EBO equations (16) have 
meaningful solution only when they satisfy the following gauge invariance condition [12] for systems 
of three electronic states, 

1 r 2n - 

— / u>(n).dn = m, m = 1,2,3, .... (17) 

27T Jo 

We choose different ratios of «(n), /3(n) and 7(11) and calculate corresponding u> s as: (a) a — f3 
= 7, u = ±Va{3 + 2sina}5; (b) a = 2/3 = 7, u = ±%{9 + 8 sin f }5; ( c ) 2a = (3 
= 27, u = ±Va{6 + 2 sin 2a} 1 : where a(n) is mixing as well as ADT angle (upto an additive 
constant) among the electronic states. It is important to note that in all such situations, divergence 
equations (12) have non - zero contributions. When a(n) is the function of two nuclear coordinates, 
x(— qcosO) and y(— gsin#)and is being equated as a(n) = |, the product, Vct(n) • dw = \dd. For 
realistic systems, a{6) can be calculated as function of 9 from the electronic eigenfunctions of the 
equation, H e (e, n)£j(e, n) = Uj(n)^(e, n). Thus, the gauge invariant integrals for the above three 
cases are, 



r\ = 


1 1 

2^io 2 1 




1 1 


r 2 = 


2^io 4 1 




1 1 


r 3 = 


2% Jo 2 



(7 f 4 

-2sin-}^0 = 4V5 / 4 (1 - -sin 2 0)5rf(/)= 1.03 (18) 
2 Jo 5 

■{9 + Ssm^de = 2Vl7 - ^sin 2 (/))5d0 = 0.934 (19) 
4 Jo 17 

t~ 1 

f2sin0}sd0 = 2^2 / 4 (1 - -sin 2 0)5^ 

jo 2 

+ 2V2 / 4 (1- -sin 2 <P)id(p= 1.216 (20) 
jo 2 

with the general form of incomplete elliptic integral of the second kind. It is quite obvious that since 
a(n) (A is analytic) is analytic, the nature of these integrals will be generic for any functional form 
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of a(n). When three electronic states are coupled, the non - adiabatic effect of the upper states on 
the ground is equivalent to a potential developed due to elliptic motion of the nuclei around the 
point of conical intersection. Moreover, single surface EBO can be derived only for specific ratios of 
mixing/ADT angles, e.g., gauge invariance condition is approximately obeyed in case (a), whereas 
in the cases (b) and (c), integrals are either away or far away from the gauge condition. 

Since the general form of f and A (equations (5) - (7)) with any arbitrary ratios of mixing angles 
satisfies the equation, VA + r A = 0, and ensures the ADT, uniquely defined diabatic potential 
matrix in configuration space is guaranteed by the unit matrix, D = exp(/ 27r r ■ dn). The explicit 
expression of D [13] is derived for the case (a) by using the corresponding G matrix (equation (15)) 
as, 



D 



r2n 

= G exp(— i / d3(n) ■ dn) 
J o 



G* 



1 

3 + 2. si 



x 



(si + l) 2 
+ (1 + cs 2 )C\ 



— (cs 2 + si) 
x(3 + 2si)^S 1 
+2cs 3 S 2 

— (cs.si — cs) 
x(3 + 2.si)^S 1 
y +2(1 + si + si.cs 2 ))S 2 

I 1 \ 
1 



— [cs 2 + si) 
x(3 + 2.si)^S 1 
+2.cs 3 S 2 

(cs.si — cs) 2 
[l + 2si + 2.si 2 
+2.si.cs 2 + cs 4 )Ci 

(si + l)(3 + 2.si)^S 1 
+2.cs(si 3 + cs 2 - si 2 )S 2 



(cs.si — cs) 
x(3 + 2.si)^S 1 
+2(1 + si + si.cs 2 )S 2 

-(si + l)(3 + 2.si)sSi 
+2.cs(si 3 + cs 2 - si 2 )S 2 



(cs 2 + si) 2 
+ (2 + 2. si 3 
+si 2 + cs 2 )Cx 



(21) 



\0 I) 

where si = sin a, cs = cos a, C\ = cos^yr]^) ~ 1, Si — sin^^rx) ~ and S 2 = sin^^rx) ~ and 
for the other two cases (b) and (c), Ci ^ 1 and Si, S 2 ^ 0. 

In equation (15), the contribution of non - adiabatic effects appear as ±iu> in the KE operator 
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(diagonal) as well as in the potential energy matrix (through G matrix). Since single surface EBO 
equation is derived by neglecting the effect of G matrix (Z)| =2 G*^{ui — u^ipj), diabatic and single 
surface EBO equations will not be equivalent irrespective of any arbitrary ratios of mixing angles. 
Alternatively, we can say that for any arbitrary ratios of mixing angles, the term, X^ =2 G*j(ui—ui)ipj, 
is not negligible and single surface EBO equations can not be derived. The condition of gauge 
invariance and the existence of uniquely define diabatic potential from non - adiabatic terms clearly 
tell that only for specific ratios of mixing angles, single surface EBO equations is acceptable. 

We justify our theory by performing numerical calculations on a two coordinate quasi -"JT scat- 
tering" model [13] in which the harmonic oscillator potential and linear coupling term are replaced 
by more general potentials. The adiabatic PESs ui, 1 = 1,2,3 are given below: 





= t^o- 


-^ l {x)) 2 y 2 + A l x f(x,y) 


u 2 (x,y) 


1 2 2 
= ^ U ° y 


- (D 1 -A 1 )xf(x,y) + D 1 


u 3 (x,y) 


1 2 2 

= 2^ u ° y 


- (D 2 -A 1 )xf(x,y) + D 2 




= o>iexp(( 






f(x,y) 


= exp( 


x 2 + y 2 ^ 
a 2 ) 



where // = 0.58 amu, A x = 3.0 eV, D x = 5.0 eV, D 2 = 10.0 eV, uj = 39.14 x 10 13 s~\ uj± = 
7.83 x 10 13 a = 0.3 A and U\ = 0.75 A, also Cartesian coordinates x and y are defined in the 
intervals — oo < x < oo and — oo < y < oo and related with polar coordinates as, x = qcos6 and 
y = q sin 9. These adiabatic potentials describe a two arrangement channel system where x — > oo 
and x — > — oo are the reagents and products asymptote, respectively. 

Introducing a(n) = | = |tan _1 (^), we construct diabatic potential matrices (W), vis -a- vis, 
diabatic and single surface EBO equations for all the three cases and then, initialize the wavefunction 
on the ground vibrational state with different initial KE at the asymptote of reagents, propagate 
the time dependent wavefunction using DVR [18] and project the final wavefunction with the 
asymptotic eigenfunctions of the Hamiltonian to calculate state - to - state vibrational transition 
probabilities. Dynamical calculations are carried out at total energy 1.20 and 1.80 eV. Since the 
point of conical intersection is at 3.0 eV, upper electronic states are expected to be classically 
closed at those energies. We demonstrate all the results in Table (I) - (III) for the cases (a) - 



9 



(c), respectively. Table(I) exhibits that reactive transition probabilities calculated by single surface 
EBO equations not only follow the correct symmetry but also achieve quantitative agreement with 
diabatic results whereas in Table (II) and (III), single surface EBO results are symmetry broken as 
well as inaccurate. Thus, these numerical results predict that only for the case (a), single surface 
EBO equation can be constructed. 

In summary, we have used the generalized form of real orthogonal electronic basis functions in terms 
of mixing angles (a(n), /3(n) and 7(n)) among the three electronic states and expressed the NAC 
terms of adiabatic nuclear SEs with these angles. When ADT angles are chosen as mixing angles, 
adiabatic nuclear SEs transform to diabatic SEs, i.e, ADT and mixing angles are equal upto an 
additive constant. Since the NAC (r), ADT (A) matrices satisfy the ADT condition and mixing 
angles are analytic, we find a curl condition is also satisfied with non - zero divergence for each 
element of f matrix. In irrotational case, solution of these curl equations imply that mixing/ADT 
angles are related with integer ratios. Consequently, both the adiabatic and diabatic nuclear SEs are 
being simplified extensively. Single surface EBO equations derived from adiabatic nuclear SEs have 
meaningful solution and quantitative agreement with corresponding diabatic case only for specific 
ratios of mixing angles where EBO equation is gauge invariant and provide uniquely defined diabatic 
potential energy matrix. We also find that the non - adiabatic effect associated with single surface 
EBO equation derived from three coupled electronic state is a potential that arises due to the elliptic 
motion of nuclei around the conical intersection. 

We acknowledge Department of Science and Technology (DST, Government of India) for 

financial support through the project no. SP/S1/H-53/01. S.A. would like to thank Professor S. 
P. Bhattacharyya and Professor J. K. Bhattacharjee, I.A.C.S., Kolkata for their comments on curl 
and divergence equations. 
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TABLE I 

Reactive state - to - state transition probabilities. Three diabatic surfaces are constructed consid- 
ering the relation, a(n) = (3(n) = 7(11) where the EBO is derived under the same situation. 



E (eV) 0^0 -> 1 0^2 0^3 0^4 0^5 0^6 

1.20 0.0229 a 0.0000 0.0616 0.0000 0.0006 

1.20 0.0279 6 0.0032 0.0656 0.0049 0.0010 

1.80 0.1002 0.0000 0.0401 0.0000 0.0912 0.0000 0.0238 

1.80 0.1194 0.0069 0.0524 0.0016 0.1196 0.0093 0.0139 



a Diabatic 
b EBO 
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TABLE II 

Reactive state - to - state transition probabilities. Three diabatic surfaces are constructed consid- 
ering the relation, a(n) = 2/3(n) = 7(11) where the EBO is derived under the same situation. 



E (eV) 0^0 -> 1 0^2 0^3 0^4 0^5 0^6 

1.20 0.0225 a 0.0000 0.0491 0.0000 0.0124 

1.20 0.0114 6 0.0321 0.0298 0.0230 0.0002 

1.80 0.1389 0.0001 0.0246 0.0000 0.0876 0.0000 0.0103 

1.80 0.0732 0.0684 0.0389 0.0036 0.0835 0.0190 0.0129 



a Diabatic 
b EBO 
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TABLE III 

Reactive state - to - state transition probabilities. Three diabatic surfaces are constructed consid- 
ering the relation, 2a (n) = /3(n) = 27(11) where the EBO is derived under the same situation. 



E (eV) 0^0 -> 1 0^2 0^3 0^4 0^5 0^6 

1.20 0.0665 a 0.0000 0.0115 0.0000 0.0008 

1.20 0.0315 6 0.0102 0.0629 0.0002 0.001 

1.80 0.1220 0.0000 0.0390 0.0000 0.0697 0.0000 0.0027 

1.80 0.1318 0.0309 0.0363 0.0064 0.0586 0.0573 0.0061 



a Diabatic 
b EBO 
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